Progressive motion of an ac-driven kink in an annular damped system 
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A novel dynamical effect is presented: systematic drift of a topological soliton in ac-driven weakly 
damped systems with periodic boundary conditions. The effect is demonstrated in detail for a 
long annular Josephson junction. Unlike earlier considered cases of the ac-driven motion of fluxons 
(kinks) , in the present case the long junction is spatially uniform. Numerical simulations reveal that 
progressive motion of the fluxon commences if the amplitude of the ac drive exceeds a threshold 
value. The direction of the motion is randomly selected by initial conditions, and a strong hysteresis 
is observed. An analytical approach to the problem is based on consideration of the interaction 
between plasma waves emitted by the fluxon under the action of the ac drive and the fluxon itself, 
after the waves complete round trip in the annular junction. The analysis predicts instability of 
the zero-average-velocity state of the fluxon interacting with its own radiation tails, provided that 
the drive's amplitude exceeds an explicitly found threshold. The result is valid if the phase shift ip 
of the radiation wave, gained after the round trip, is such that sin^ < 0, the threshold amplitude 
strongly depending on (p. A very similar dependence is found in the simulations, testifying to the 
relevance of the analytical consideration. 

PACS numbers: 05.45.Yv 85.25.Cp, 74.50. +r, 
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I. INTRODUCTION 

It is well known that progressive motion of topological 
solitons (kinks ) in lossy media can be supported by a 
dc (constant) field which couples to the soliton's topo- 
logical charge Q. In this case, the kink is, effectively, 
a quasiparticle. The same quasiparticle approximation 
strongly suggests that ac (time-periodic with zero mean 
value) field(s) applied to the kink in an infinitely long 
homogeneous system can only give rise to its oscillatory 
motion. 

A more sophisticated possibility is to drive a kink by 
a pure ac force in a lossy system subjected to a periodic 
spatial modulation. In particular, spatial periodicity is 
an inherent feature of discrete systems (dynamical lat- 
tices). As it was first predicted in an analytical form 
for damped lattices of the Frcnkcl-Kontorova (FK) and 
Toda types in Refs. and ||, respectively, and for 
a periodically modulated lossy continuum system of the 
sine-Gordon (sG) type in Ref. Q, resonances (of differ- 
ent orders) between the time-periodic driving force and 
periodic passage of a moving kink through the spatial 
inhomogeneities make it possible to support progressive 
motion of the kink without any dc field. This mecha- 
nism selects absolute values of the mean velocity of the 



ac-driven kink which provide for the resonance of the or- 
der m : n, 
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where u) and L are the frequency of the ac drive and 
the period of the spatial modulation. The sign of the 
velocity is determined by a random initial push applied 
to the kink. While the velocity selected by the locking 
condition ([!]) does not depend on the amplitude of the 
ac drive, the driven progressive motion cannot take place 
unless the drive's amplitude e exceeds a certain minimum 
(threshold) value ethn which is usually proportional to 
the friction coefficient ||, U . 

The effect was generalized to include a case when a 
combination of ac and dc fields is applied [||. In this 
case, a constant-velocity step is predicted, in the form of 
an ac-frequency-locked value (]l|) of the mean velocity of 
the driven kink in a finite interval of values of the dc 
component of the drive. 

The possibility of stable ac-driven motion of the soliton 
in the lossy Toda lattice has been later demonstrated in 
simulations Q and in a laboratory experiment with an 
electrical transmission line described by this model [fil. 
Search for such regimes in FK models turned out to be 
more difficult, but simulations have finally demonstrated 
that the soliton in the form of a dislocation may be driven 
with a nonzero mean velocity by the ac force || . 

In terms of physical applications, the most appropriate 
medium where these effects may be realized is provided 
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by long Josephson junctions (LJJs) j|, ||, which sup- 
port topological solitons in the form of trapped fluxons 
(magnetic- flux quanta ||). A convenient way to imple- 
ment the periodic spatial modulation in LJJ is to apply 
constant (dc) magnetic field to an annular (and suffi- 
ciently long) junction, with the field's vector lying in the 
plane of the annular junction. In this case, the full length 
of the ring-shaped junction plays the role of the modu- 
lation period L in Eq. (Q). Note that the use of the LJJ 
of the annular shape is very natural also because it pre- 
serves the number of fluxons trapped in it. The necessary 
drive can be readily provided by ac electric bias current 
distributed along the junction. 



As it was first proposed in Ref. |10|, this configuration 
gives rise to a harmonic spatial modulation in the cor- 
responding sG model, the modulation amplitude being 
proportional to the external magnetic field (the pertur- 
bation generated by the magnetic field is formally tan- 
tamount to an additional spatially modulated dc drive 
with zero average). The resonant relation (Q) applies to 
this case, the period L being the circumference of the 
annular junction, as it was mentioned above. Experi- 
mental observation of the ac-driven motion of a fluxon in 
the annular LJJ in the presence of the constant magnetic 
field was reported in Ref. , in the form of nonzero dc 
voltage across the junction, induced by ac bias current 
applied to it. 

A related issue is progressive motion of the fluxon un- 
der the action of a pure ac drive, or even a random drive, 
in the so-called ratchet (asymmetric) effective potential 
|Q. However, a principal difference is that, in the case of 
the ratchet, there is an a priori selected preferred direc- 
tion of motion. Experimental realization of a ratchet for 
fluxons in LJJs was very recently reported in Ref. [fL3| . 

Thus, ac-driven progressive motion of a topological 
soliton in an indefinitely long lossy system can only be 
possible if it is assisted by the spatial periodic modula- 
tion. The objective of this work is to show, both nu- 
merically and analytically, that, quite surprisingly, the 
ordinary spatially uniform ac drive may support progres- 
sive, on average, motion of a fluxon in either direction 
in the perfectly uniform annular junction with non-zero 
losses, provided that the amplitude of the drive exceeds 
a certain threshold value. The effect is clearly stipu- 
lated by the finite size of the ring-like junction, being 
possible because a fluxon perturbed by the driving field 
may emit small-amplitude (quasi-linear) waves (radia- 
tion, alias "plasmons" or "phonons") propagating along 
the junction and eventually hitting the same fluxon. 

We note that a somewhat similar effect was observed 
in numerical simulations reported in recent preprints 
Jl4| , [T^ |, which, however, was related to the above- 
mentioned ratchet mechanism supporting unidirectional 
motion of ac-driven fluxons, and the ac drive was always a 
two- frequency one. It was stated that the progressive mo- 
tion is possible due to the simultaneous action of driving 
forces oscillating at frequencies u> and 2u> in the absence 
of static spatial potential. We stress that in our work 



the ac drive is always represented by a single-frequency 
harmonic. 

The rest of the paper is organized as follows. Numer- 
ical results which show the ac-driven progressive motion 
of the fluxon are displayed in section H. A theoretical ex- 
planation, that reflects our understanding of the effect, 
which is based on interaction of the fluxon with its own 
radiation "tails" , is set forth in section III . Comparison 



of the analytical and numerical results is presented in 
section IV , and section [v] concludes the work. 



II. NUMERICAL RESULTS 

A. The model and computational procedure 

The ac-driven weakly damped annular Josephson junc- 
tion is described by the well-known p_erturbed sG equa- 
tion for the Josephson phase 4>(x, t) 



sin 



4> = a(j} t - 7a C sin(a;i), 



(2) 



where the subscripts stand for the partial derivatives, 
the coordinate x, running along the junction, and the 
time t are normalized so that the Josephson penetration 
length and plasma frequency are both equal to 1, a is the 
damping coefficient, and 7 ac is the amplitude of ac drive. 
In the infinite system without perturbations (a — 7 ac = 
0), the fluxon is described by the kink solution to the sG 
equation, 



ifi = 4 arctan 



exp a 



(3) 



where u is the velocity of the fluxon, a — ±1 is its polar- 
ity, and £ is the instantaneous coordinate of its center, 
which is = £ + ut in the case of free motion of the 
fluxon. If a single fluxon is trapped in the annular junc- 
tion, Eq. (^) is supplemented by the boundary conditions 
(b.c.) 



<KL,t) 
t>x(L,t) 



0(0, t) + 2tt, 
0*(O,i), 



(4) 
(5) 



where L is the circumference of the ring. 



Equation (|) with b.c. (g) and (§) was simulated by 
means of the "StkJJ" software package, described in de- 
tail elsewhere jl6|. The results have been checked by 
comparing them with simulations carried out with the 
"Soliton" package jl7|, and good agreement was found 
in all the cases considered. 

The simulations were performed as follows. We started 
with an initial state corresponding to a fluxon trapped 
in the system at the zero amplitude of ac drive. Then 
the amplitude was increased by small steps A7 ac , each 
time calculating the average dc voltage Va c across the 
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junction, 

Vdc = J dt J dx(f> t (x,t) 

= -^if dx l<t>(x,T)-<t>(x,0)], (6) 

where T is a sufficiently large time-averaging interval [in 
physical units, the average voltage is given by the expres- 
sion (|^) multiplied by the magnetic ffux quantum <3?o]- 
The value of T was taken as multiple of the ac-drive's 
period 2tt/lu. 

The dependence Vd c (7ac) characterizes the average ve- 
locity of the progressive motion of the fluxon in the an- 
nular junction. Indeed, substituting a solution for the 
moving fluxon in the form <j}{x,t) — 4>n(x — £(t\), see 
Eq. (||), into the first integral expression in Eq. (g), we 
obtain 

- 

where b.c. (||) was used. A natural definition of the 
nuxon's average velocity being 

u= tJ dt d-t> 

Eq. (0) yields an expression for the average velocity in 
terms of the average dc voltage, 

u = LV dc . (8) 

The results of the simulations are displayed below in the 
form u = u(7ac), using numerically found dependencies 
Vdc(7ac) and the relation (||). This way of the presenta- 
tion of results is the most appropriate one if the objective 
is to look for dynamical regimes with a nonzero average 
velocity of the ac-driven kink. 

For some values of 7 ac the motion of fluxon is appar- 
ently chaotic (see below). In such cases, the average ve- 
locity calculated with Eq. ^ will not converge as T — > oo. 
To avoid this problem, we used two different averaging 
procedures: 

• (i) The averaging was made simply over one period 
2n/ui of the ac drive, and the amplitude 7 ac was 
changed in extremely small steps, Sja, c = 0.0001, 
while, otherwise, Sj ac = 0.01 is used. This, in a 
sense, means that instead of one point on the u(7 a c) 
dependence we get 100 closely situated points. If 
the dynamics is chaotic those 100 points will have 
large spread in u. On the other hand, all the points 
will give almost equal values of u if the dynamics is 



regular and averaging converges. Although some 
artifacts may appear around bifurcation points, 
this technique provides very good insight into dy- 
namical states of the system. 

• (ii) The averaging was performed iteratively for 
gradually increasing value of T until the conver- 
gence of u with the accuracy 0.001 was reached, 
but not longer than for T = 1000. In this case, 
for the non-chaotic states the averaging gives the 
true value of equilibrium average velocity u, while 
for chaotic states averaging was interrupted at T = 
1000, and the corresponding non-converged value 
of u was included into the plot. 

Results produced by means of both methods will be 
displayed below. 

B. The ac-driven progressive motion 

First, we present dependencies w(7ac) generated by the 
averaging method (i) . Four characteristic examples of the 
dependencies, obtained by increasing (gray points) and 
decreasing (black points) the drive's amplitude 7 ac for 
the annular LJJ of the length L — 10, with the dissi- 
pation constant a = 0.1 and four different values of the 
driving frequency lo, are shown in Fig. [j]. In Fig. 0(a), 
corresponding to ui = 0.75, the average velocity of the 
fluxon remains zero in the region 7 ac < 0.38. However, 
in the interval 0.38 < 7 ac < 0.6 the average velocity is 
different from zero, which means that the fluxon moves, 
on average, in one direction. As one sees in Fig. [I] (a), the 
signs of the average velocity are different for the branches 
of the characteristic u(7 ac ) corresponding to "sweeping 
up" and "sweeping down", i.e., increasing and decreas- 
ing 7 ac . Due to the symmetry of the system, these signs 
are actually picked up randomly, depending on the ini- 
tial conditions and peculiarities of the simulations, such 
as the way 7 ac was varied. 

At the frequency uj = 1.05, corresponding to Fig. |l|(b), 
the fluxon is set into motion in the region 7 ac > 0.1. The 
average velocity u increases with 7 ac , but at 7 ac ss 0.26 
a transition to a chaotic-like behavior takes place. In 
the latter case, the averaging over a long time interval 
does not produce any definite value of the mean velocity. 
Closer inspection of this region shows that we are deal- 
ing not with true chaos, but rather with a deterministic 
quasi-periodic regime with a nonzero average velocity. In 
particular, the fundamental frequency of the correspond- 
ing time series is, effectively, incommensurate with the 
ac driving frequency in this case, that is why the average 
velocity cannot be effectively calculated using the algo- 
rithm described above. 

Continuing the consideration of Fig. 0(b), we notice a 
split branch in the interval 0.48 < 7 ac < 0.54, where the 
average velocity is u « 0.3. Accurate inspection shows 
that the averaging during one period yields a point be- 
longing to one branch, while averaging over the next pe- 



4 



! -5 -| 
.0 
i.5 
i.O 

! .5 -| 
.0 
i.5 




! -5 -| 
.0 
! -5 -I 



(a) w=0.75 



1 *tf:-s£ I I *V," 



— 1 1 — 

(c) (o=1.12 




(d) o)=1.20 



— 1 1 1 1 1 1 1 r— 

0.0 0.2 0.4 0.6 0.8 

amplitude of ac bias current, / 



1.0 



FIG. 1: The dependences w(7ac) obtained by sweeping 7 ac 
up (black points) and down (gray points) at four different 
frequencies of the ac drive: ui — 0.75 (a), u> = 1.05 (b), ui — 
1.12 (c), and u — 1.20 (d). The length of the annular junction 
is L — 10, and the damping coefficient is a = 0.1. 



riod produces a point situated on the other branch, and 
so on. This means that the fundamental frequency of 
the oscillatory component of the fluxon's motion is half 
the driving frequency. While generation of higher har- 
monics is common to all nonlinear systems, the appear- 
ance of sub- harmonics is a more specific effect, implying 
the existence of a certain form of parametric instabil- 



ity. For larger values of 7 ac , the system enters the region 
of chaos, which is shortly interrupted by a deterministic 
branch with a nonzero average velocity in the interval 
0.75 < 7 ac < 0.80. 

For driving frequencies that are ~ 10% — 20% higher 
than the Josephson plasma frequency (which is 1 in the 
notation adopted above), one can observe the fluxon's 
motion with a nonzero average velocity in a broad range 
of 7ac, see Figs. [j](c) and (d). It is noteworthy that the 
average velocity may attain relatively large values, which 
become comparable to the Swihart velocity of the junc- 
tion (= 1 in our notation). For instance, in Fig. |l|(d) 
u w 0.5 at 7 ac = 0.5, and reaches the value w 0.7 for 
larger values of 7 ac . 

It is also interesting to note that the dependence 
w(7ac) shown in Fig. |l](d) features an additional step at 
7ac ~ 0.6, which suggests that there are two different 
mechanisms which generate the nonzero average veloc- 
ity. The additional step may correspond to switching of 
the system between these two different mechanisms. 

To visualize the fluxon motion, in Fig. we present 
the law of motion, £(t), of the fluxon's center for differ- 
ent values of the ac driving amplitude and fixed driving 
frequency u> = 1.05, which corresponds to the dependence 
u(7ac) shown in Fig. 0(b). It is clearly seen in Fig. ||(a) 
that, indeed, at 7 ac = 0.2 a systematic drift of the fluxon 
in one direction is superimposed on periodic oscillations. 
In Fig. H (b), which corresponds to the quasi-periodic 
motion at 7 ac = 0.4, we see that the fluxon motion is 
qualitatively similar to that displayed in Fig. 0, but it is 
less regular and nonperiodic. 

The fluxon's law of motion for the case which corre- 
sponds to the generation of the sub-harmonic at 7 ac = 0.5 
is shown in Fig. 0(c). One can easily see that the period 
of the oscillatory part of the soliton's motion in this case 
is indeed twice the period of the ac drive, as the trajec- 
tories look slightly different during alternating (even and 
the odd) periods of the ac drive. 

Finally, in Fig. ||(d) one can see the trajectories cor- 
responding to the point 7 ac = 0.8, which belongs to a 
narrow window in Fig. 0(b), where a deterministic law 
of motion becomes stable again. Qualitatively, this plot 
is similar to that in Fig. f|(a), but has a higher average 
velocity. 

Results produced by the averaging method (ii) are dis- 
played, for to = 1.05, in Fig. ||(a). As is seen from the 
comparison of this figure with Fig. 0(b) , details of the pic- 
tures generated by the averaging procedures (i) and (ii) 
are somewhat different inside the chaotic regions. Nev- 
ertheless, both methods identify the same regions of the 
regular motion, and yield the same dependencies 7t(7 ac ) 
in the regular cases. In fact, the consideration of the reg- 
ular motion with nonzero average velocity, rather than of 
chaotic regimes, is the main subject of this work. 

Our numerical simulations of the perturbed sG equa- 
tion (0) used a finite-difference scheme. Since the ac- 
driven progressive motion of kinks is possible in indefi- 
nitely long discrete lossy lattices |2| 0, ^, |^| , it is necessary 
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FIG. 2: The coordinate of the fluxon's center as a function of time for the ac-drive's frequency u = 1.05 [which corresponds 
to the case shown in Fig. |l](b)] and different values of the drive's amplitude 7 ac : 0.2 (a), 0.4 (b), 0.5 (c), and 0.8 (d). The 
rightmost plot shows the driving ac signal sin(ui) vs. time. 



to check that the effects reported above are not artifacts 
produced by the discretization of the model. To this end, 
simulations were run with the smaller value of the step 
size Ax of the finite-difference scheme. Comparing the 
results, we have concluded that they are almost iden- 
tical, except for small variations in the location of the 
thresholds at which the system switches from one state 
to the other. For example, Fig. ||(b) displays the depen- 
dence w(7ac) obtained for the same case as in Fig. |l|(b), 
but with Ax = 0.005, rather than Ax = 0.01 used in 
Fig. [l](b) . As for the different signs of the average veloc- 
ity in some regions of the regular motion in Figs. 0(b) 
and |(b), this sign is always randomly selected by the 
initial conditions, as it was explained above. 

The numerical results presented above clearly demon- 
strate that progressive motion of fiuxon can indeed be 
supported, in either direction, by the uniformly dis- 
tributed ac drive (bias current) in annular uniform 
weakly damped LJJ. An explanation to the effect will 
be proposed in the next section. 

III. ANALYTICAL CONSIDERATION 

A. Radiation waves emitted by the ac-driven 
fiuxon 

The analysis developed below is based on a basic prop- 
erty of the system under consideration: the fiuxon oscil- 
lating under the action of the ac drive emits plasma waves 
in both directions, and, due to the annular shape of the 
junction, these waves interact with the same fiuxon af- 
ter completing a round trip. However, this property does 
not offer an immediate explanation to the possibility that 
the fiuxon will be able to systematically drift in one di- 



rection, if pushed initially. Therefore, a detailed analysis 
is necessary, which is developed below. 

As it is known from earlier works on the perturbation 
theory (PT) for the sG equation (see a review 0]) , a kink 
(fiuxon) oscillating, without any systematic progressive 
motion, under the action of the ac drive emits radiation, 
in the lowest approximation, to the left and to the right 
at two wavenumbers ±fc, where 

k = y^ 2 - 1, (9) 

provided that u> > 1, i.e., that the driving frequency ex- 
ceeds the plasma frequency of the junction. If lj < 1, the 
emission takes place in higher orders of PT (at higher 
harmonics). 

Here, we concentrate on the case u> > 1. The numeri- 
cal results (both displayed in Fig. [j] and those which are 
not shown here) suggest that the strongest effect takes 
place when (u 2 — l) is a positive but relatively small pa- 
rameter. In this case, the application of PT, which needs 
the drive's amplitude 7 ac and the damping coefficient a 
to be small, predicts that the oscillating fiuxon emits 
two waves, which, far from the fiuxon, have the following 
form (neglecting, for the time being, the action of the 
dissipation): 

<p± = A a sin (±kx - u)t) . (10) 

The amplitude of the two emitted waves is predicted by 
PT to be § 

A a = 7 ac/fc 2 , (11) 

which is correct if 

k 2 > 7ac /3 and fc2 > ( 12 ) 
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FIG. 3: The dependence u(y a c) produced using: (a) the al- 
ternative averaging procedure, cf. Fig. |l|(b) ; and (b) a smaller 
spatial step of the finite-difference integration scheme. The 
scattered points in panel (a) are those at which the averag- 
ing did not converge, thus they mark regions of the chaotic 
motion of the ac-driven fluxon. 



In fact, the former inequality does not hold in many cases 
when the simulation reveal the progressive motion of the 
ac-driven fluxon, but this is not a principal limitation: 
instead of the expression (|TT|) , one can use more general 
ones (6.23) and (6.18) from Rcf. jlj]. In such a case, 
a relation between Aq and the drive's amplitude 7 ac is 
more complex than that given by Eq. (|ll] ) , but there will 
be no essential change in the analysis presented below, 
nor significant changes in the final results. 

Attenuation of the emitted waves due to the presence 
of the loss term in Eq. (j^) will play a crucially important 
role below. To consider the attenuation, it is convenient 
to use a reference frame (x',t') moving at a constant 
velocity u, which, as well as in the previous section, is the 
average velocity of the fluxon (if it moves on average). In 
the analysis, we will treat u as the smallest parameter, 
in comparison with 7 ac , a, and lo 2 — 1 = fc 2 , the main 
objective being to show that, at a certain threshold value 
(7ac)thr 01 the ac-drive's amplitude, the state with u = 
becomes unstable, and a symmetric pair of new nontrivial 
states with finite but very small average velocities ±u 
appears as a result of an instability-triggered pitchfork 
bifurcation if 7 ac slightly exceeds 7th- 

In the moving reference frame, the linearized version of 
Eq. (13), which governs the propagation and attenuation 



of the radiation waves, becomes 



0: 



(13) 



where x' = (x — ut) / yl — u 2 and t' = (t— ux) \J\ — u 2 . 
A stationary shape of the radiation wave in the moving 
reference frame is sought for as [cf. Eq. (flO|)] 



i>± = A ± (x') sin (±\/ 



l-x 1 



(14) 



where a slow dependence A{x') is produced by the dis- 
sipative damping of the wave, and lo' is the driving fre- 
quency in the moving reference frame, 



(15) 



In the first approximation, which assumes that both the 
loss constant a and the average velocity u are small, the 
substitution of Eq. ( |l4| ) into Eq. (|l3|) gives rise to an 
equation 



dA ± 
"dx 7 



(16) 



where the notation for the radiation-wave's group veloc- 
ity was introduced 

did' 



dk7 = VuJ^l/cj'^k[l + u 2 /(2k 2 )]. (17) 

To obtain the final expression in Eq. ( |l7j ) , it was assumed 
that, in accord with what was said above, « 2 < fc 2 < 1 
(hence the second term in the square brackets is a small 
correction) . 

A solution to Eq. (|q) is obvious: 



A±(x') = A exp 



2u'^y/ 1 - u 2 



(l±u' SI u)\x'\ 



(18) 

In particular, we will need the value of the emitted-wave's 
amplitude after the completion of the round trip, from 
x' = to x' = L' = Lyl — u 2 , where the Lorentz con- 
traction is taken into regard. According to Eqs. Jlg| ) and 
(0), we find 



A±{L')=A ( , 



(19) 



where 

A = A exp 



aL 
2uL 



4fc 3 



(20) 

(hereafter, it is implied that k = \/ui 2 — 1 is taken as a 
positive square root). As it was said above, the funda- 
mental assumption in the analysis is that the velocity u 
is very small. Nevertheless, small corrections ~ u and 
~ u 2 , which are retained in Eqs. ( JToj ) and (20), respec- 
tively, will play an important role below. 
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B. Dragging the fluxon by the radiation waves in 
the circular system 

As we are dealing with the circular system, each emit- 
ted wave, after having completed its round trip, strikes 
the fluxon and exerts some dragging force acting on it. 
Dragging a kink by a radiation wave passing through it 
was analyzed in detail in Ref . p9|| . It was found that a 
nonzero dragging force appears at the second order in the 
dragging-wave's amplitude, provided that a^O. An ex- 
pression for the dragging forces induced by the two waves 
(j>± can be shown to take the following form: 



It is straightforward to derive in the lowest (adiabatic) 
approximation jj] an equation of motion for the coordi- 
nate £ of the fluxon driven by the terms on the right-hand 
side of Eq. (pi|). In the "non-relativistic' approximation 
(for small velocities), it takes the form 



(7ac/fc 2 ) sin (ut) 



±(tt/4) 2 akA 2 ± 



(21) 



(this expression assumes that the kink's mass is normal- 
ized to be 1). However, taking the expressions © for the 
amp litudes A± and substituting them into the expression 
(pl|), one can readily see that the resultant force, pro- 
duced by the asymmetry between the two waves, would 
be not accelerating, but braking the fluxon's motion. 

A key point for the explanation of the possibility of the 
ac-driven progressive motion is to notice that the waves 
( |l4| ) act on the fluxon in combination with the direct 
ac drive. To take all the forces into regard, a known 
perturbative technique may be applied |l[ : the net wave 
field is sought for in the form 



d 2 £ d£ _ (77T 

~d¥ + a !t ~ T 

+A [sin (ut -kL- kg + sin (ut - kL + kg)] 
——aLA u [sin (ut — kL — kg) — sin (ut — kL + kg$j) 

[recall that a = ±1 is the fluxon's polarity, see Eq. (|J)]. 

The right-hand side (r.h.s.) of Eq. ( p5| ) may be ex- 
panded for small g Then one concludes that the first 
term and the following pair of terms on r.h.s. are essen- 
tially similar to each other (in the lowest approximation, 
they do not contain g), but, with regard to Eqs. ( JTT| ) and 
(g0|), the amplitude in front of the latter pair of terms, 
Aq, differs from the amplitude in front of the first term 
by the exponential factor exp [— <xLj (2k)], that we as- 
sume to be small enough. Therefore, this pair may be 
neglected, and, keeping a term linear in £ which is pro- 
duced by the expansion of the last pair of terms on r.h.s. 
in Eq. (p5|), we arrive at a simplified equation of motion 
for the driven fluxon, 



(x,t) = ^(x,t)+<l> (} (t)+<l> + (x,t)+ ( j>_(x,t), (22) ^| +a ^i = ^ 

dt 2 



dt 



an 
T 



-^p sin (ut) + aLAoku cos (ut — kL) £ 
k l 



where tMx, t) is the field of the fluxon, <f>± are the emitted 
waves (|l(]|), and 



)(*) 



huc/k 2 ) sm(ut) 



(23) 



(26) 

We seek for a solution to Eq. ( P6|) in a natural form 
(similar to that employed in Ref. p9[ ), 



is a term representing uniform background oscillations 
generated by the ac drive, as it follows from a straightfor- 
ward solution to the linearized equation (|^) (recall that 
k = \/u 2 — 1). Then, an effective equation for the fluxon 
field ib is obtained by the substitution of the combina- 
tion ( p2] ) into the underlying perturbed sG equation (||), 
expansion of sin <f>, and taking into regard the expression 
( p3| ) for the background component of the field. The final 
result (written in the laboratory reference frame) is 

iptt - ipxx + sin ip + atpt = - (1 - cos ip) • { 
Aq [sin (kx — ut + kL) + sin (—kx — ut + kL)] 

-aLAoii [sin (kx — ut + kL) — sin (—kx — ut + kL)] 

(24) 



4fc 2 



■ sin (ut) + fo(*)) 



(27) 



1 



2 

(7ac/fc 2 ) sin (ut)}, 



where the amplitude Aq is defined by Eq. (|2(]|), and the 
corrections =f( 1/2) aLu from Eq. ( |l9| ) are taken into ac- 
count. The phase shift kL appearing on the right-hand 
side of Eq. ( |2~4| ) is accumulated by the wave after the 
completion of the round trip along the ring. Note that 
the renormalization of the driving amplitude in the first 
term on the right-hand side of Eq. (^4]) is quite essential 
in the present case, when k 2 = u 2 



where the first term is a response to the ac driving force 
on r.h.s. of Eq. ( p6[ ) (we neglect a small correction 
to it from the friction force on the left-hand side of the 
equation, and set 1/u 2 ~ 1, which is true in the case 
under consideration), and the term £o(t) takes into regard 
a possibility of a slow systematic motion (drift) of the 
fluxon. 

In the next approximation, we replace the multiplier 
in the second term on r.h.s. of Eq. ( p^ ) by the first 
term from the expression (p7j). In order to single out 
terms contributing to the slow drift, we perform averag- 
ing in the rapid oscillations. Thus we arrive at an effec- 
tive evolution equation for the slow variable £o> in which 
<i£o /dt may be identified by the average velocity u of the 
systematic motion of the fluxon. Finally, the equation for 
the slow motion can be cast into the form of a first-order 
evolution equation for the average velocity, 



dll _ 7T 2 LA ^ ,, T . 

-au tttt, (cm) sm [kL) 



dt 



32k 



(28) 



1 is small. 



The amplitude A in Eq. (^8|) can be replaced by the 
expressions ( |TT| ) and (^Oj), which casts the equation into 
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a form 

du 
~db 



(7ac)' 
7'hr 



au 



1 / irLa^ 

2 V 8fc 3 



exp 



~ -in (/•■£) 12!) I 



where the threshold value of the drive's amplitude is de- 
fined as 



7thr = 



(32/tt 2 ) k 3 (ah 
L sin (kL) 6XP V 2fc 



(30) 



The equation of motion ( |29| ) has a trivial stationary 
solution (trivial fixed point, FP) u = 0, which, obviously, 
corresponds to no systematic motion of the fluxon. A 
crucially important issue is the stability of this FP. The 
linearization of Eq. (29) for u — > immediately demon- 



strates that the trivial FP is unstable if the drive's am- 



plitude exceeds the threshold value given by Eq. (30). 
In this case, a transition to a nontrivial regime with 
a nonzero velocity of the progressive motion (in other 
words, a bifurcation) must take place. 

An essential feature of the expression (|30|) is its depen- 
dence on the phase kL gained by the emitted wave as a 
result of the round trip. In fact, the ac-driven motion is 
predicted to exist only if 7 t 2 hr > 0, or, as it follows from 
Eq. ©, 



sin (kL) < 0. 



(31) 



In the opposite case, the effect is absent, or, possibly, the 
threshold amplitude is very large (then the perturbative 
approach is irrelevant). 

If 7 2 c > 7j hl . > 0, the existence of two mutually sym- 
metric nontrivial FPs, with finite velocities ± u, may be 
expected. One can try to find these velocities close to the 
bifurcation point, i.e., for a case < 7 2 c — 7 t 2 hr <C 7 t 2 hr , 
taking into account the cubic corrections in Eq. (p^). In 
fact, in the particular approximation in which the cubic 
term, which is present in Eq. (p9[), was derived, a formal 
result will be u 2 < 0, which, actually, implies that the 
bifurcation is subcritical, i.e., it takes place by a finite 
jump; in the opposite case, one should expect a soft su- 
percritical transition (corresponding to a bifurcation of 
the usual pitchfork type p8[) without a jump. 

The bifurcations observed in Figs. 1(a) and, especially, 
in Fig. |l|(b) look very much like a supercritical pitchfork 
bifurcation, while Figs. 1(c) and 1(d), corresponding to 



larger values of the parameter k 



1, strongly sug- 



gest that a subcritical bifurcation (the one giving rise to 
a jump) takes place in these cases. In many other runs of 
the simulations, not displayed here, we observed a very 
similar trend, namely, to have a supercritical bifurcation 
for very small values of k 2 , and a transition to a sub- 
critical bifurcation at somewhat larger k 2 . An accurate 
prediction of the type of the bifurcation (sub- or super- 
critical) is a rather difficult issue, as it demands precise 



summation of the lowest-order nonlinear corrections at 
the order ?J 2 , that may originate from many different 
terms in the above analysis. 

Lastly, we note that the case shown in Fig. |l|(a), with 
the driving frequency ui = 0.75, which is smaller than 1, 
cannot be directly explained in the framework of the ap- 
proach developed here. However, it seems very plausible 
that the persistent motion of the ac-driven fluxon in this 
case can be explained if one takes into regard the fact 
that the kink oscillating under the action of the driving 
frequency ui < 1 emits radiation at the second-harmonic 
frequency 2u>, provided that 2a; > 1 jjj, which is the 
case in Fig. |l|(a). In accord with this, the effect seen in 
Fig. [|(a) is very weak because it is accounted for by the 
second order of PT. 



IV. COMPARISON OF ANALYTICAL AND 
NUMERICAL RESULTS 

The main prediction of the analysis presented in the 
previous section is the loss of stability of the zero-average- 
velocity state of the fluxon in the ac-driven weakly 
damped annular uniform LJJ. As a result, the fluxon 
must inevitably perform a transition to a nontrivial state 
with a nonzero average velocity, which we consider as the 
cause for the effect revealed by the simulations. Besides 
that, the analysis may, in principle, predict the character 
of the bifurcation - super- or sub-critical. 

While the comparison of the particular type of the ana- 
lytically predicted bifurcation with results of the simula- 
tions is a sophisticated problem, it is much more straight- 
forward (and more important) to compare the theoreti- 
cally predicted and numerically found points of the tran- 
sition to the progressive motion, i.e., as a matter of fact, 
dependencies of the threshold amplitude of the ac drive 
on various parameters of the system. A crucially impor- 
tant prediction of the analysis is the fact that the transi- 
tion to the nonzero velocity occurs (within the framework 
of the perturbation theory) only in those intervals of val- 
ues of the ring's length L (if all the other parameters are 
fixed) where the condition (^) holds. Another impor- 
tant and, in the same time, simple prediction is that the 
threshold value ( j30| ) remains finite (neither vanishes nor 
diverges) in the limit a — > 0. 

Comparing the analytical and numerical results, we fo- 
cused on these two basic features. In Fig. ||, the thresh- 
old value of the ac-drive's amplitude, as found from the 
simulations, is plotted vs. the length L of the annular 
junction for fixed values u) = 1.12 and a = 0.1. Salient 
peculiarities of the dependence are existence of two min- 
ima at 



\ / num \ / ] 



20.5. 



(32) 



On the other hand, the value of the wavenumber cor- 
responding to u> — 1.12 is k sa 0.53 [see Eq. ([}])], and 
Eq. (pfJ) then predicts minimum values of the threshold 
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FIG. 4: The minimum (threshold) value of the amplitude 
of the ac driving field, at which the progressive motion of 
the fluxon commences, vs. the length of the annular Joseph- 
son junction, as found from the numerical simulations. The 
driving frequency and dissipative constant are ui — 1.12 and 
a = 0.1. 

at the points where sin (kL) = — 1, i.e., at 

(33) 

Comparison of the numerical results (|32|) with the analyt- 
ical predictions ( J33| ) demonstrates very good agreement. 

Figure ^ also shows that 7 t h r diverges at some point 
between the two minima, and in a part of the region 
between them, where the data is absent in the figure, 
the effect does not take place at all (or the threshold 
is extremely high). This feature also agrees well with 
the theoretical prediction (|3(]), which yields divergences 
at points L = irn/k with integer n, and nonexistence 
of the effect in the intervals where the condition (|3ll ) 
does not hold. However, unlike the minima points, de- 
tailed comparison of the numerically found and analyti- 
cally predicted positions of the divergence points is not 
relevant, as the perturbation theory which was used to 
derive Eq. (^0|) clearly breaks down in a vicinity of the 
divergence points. 

To verify the analytical prediction for the threshold 
amplitude of the ac drive as a function of the dissipative 
constant a, the dependence 7thr(a) is shown in Fig. ||. 
The first noteworthy feature of the plot is strong hystere- 
sis. For the comparison with the analytical prediction, 
the relevant branch is the upper one, which corresponds 
to "sweeping up" , i.e., gradual increase of 7 ac , as precisely 
in this case we can detect the point at which the zero- 
average-velocity state becomes unstable, and the fluxon 
starts its progressive motion. 

As it is evident in Fig. 0, this branch of the dependence 
of 7thr vs. a indeed yields a finite value in the limit 
a — » 0, as it is predicted by Eq. (|3C|); it should be stressed, 
though, that accurate simulations are quite difficult for 



FIG. 5: The threshold value of the amplitude of the ac 
driving field vs. the dissipative constant a. The length of the 
annular Josephson junction and driving frequency are L = 10 
and ui = 1.12. The upper branch of the plot corresponds to 
the case when simulations start with a quiescent fluxon, and 
the drive's amplitude is gradually increased till the fluxon 
starts to drift. The lower branch corresponds to the opposite 
case, when the simulations start with a fluxon moving under 
the action of a sufficiently strong drive, and the amplitude is 
gradually decreased till the fluxon ceases to drift. 



very small a, as relaxation of the dynamical regime to its 
established form is very slow in this case. A particular 
value that Eq. ( |30| ) yields for a = in the case shown in 
Figs. | and [| i.e., L = 10 and lj = 1.12, is 7thr (a = 0) « 
0.2, which roughly agrees with what can be read off from 
Fig. |. 



V. CONCLUSION 

This work presents a novel dynamical effect: system- 
atic drift of a topological soliton in an ac-driven weakly 
damped system with periodic boundary conditions. The 
effect has been demonstrated in a physically relevant 
model of a long Josephson junction in the form of a ring, 
where fluxons play the role of kinks, and the ac drive is 
realized as bias current uniformly applied to the junc- 
tion. Unlike earlier considered cases of progressive mo- 
tion of the ac driven fluxon, in the present case the long 
junction and the ac driving force are spatially uniform. 
Numerical simulations reveal that progressive (on aver- 
age) motion of the fluxon commences if the amplitude 
of the ac drive exceeds a certain threshold value. With 
the increase of the amplitude, both regular and chaotic 
dynamical regimes are observed. The direction of the 
progressive motion is randomly selected by initial condi- 
tions, and regular dynamical regimes are characterized 
by strong hysteresis. The simulations demonstrate the 
effect in a well-pronounced form in the case when the 
driving frequency exceeds, but its rather close to, the 
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plasma frequency of the junction. 

The analytical approach to the problem was based on 
consideration of the interaction between plasma waves 
emitted by the fluxon oscillating under the action of the 
ac drive and the fluxon itself, after the waves complete 
the round trip along the annular junction and hit the 
fluxon. In particular, a weaker effect, which is observed 
in the case when the driving frequency is smaller than the 
plasma frequency, may be explained by the emission of 
the plasma waves at the second harmonic by the oscillat- 
ing fluxon. The main finding which the analytical consid- 
eration yields is possible instability of the zero-average- 
velocity state of the fluxon interacting with its own ra- 
diation tails. The instability sets in if the drive's ampli- 
tude exceeds an explicitly found threshold. The analysis 
predicts that the effect is only possible if the phase shift 
if = kL of the radiation wave, gained after the round trip, 
is such that sin ip < [see Eq. (pnf)], and the threshold 
amplitude strongly depends on ip, see Eq. (|3^). Numer- 
ical results show a similar dependence, and analytically 



predicted values of the length of the annular junction 
at which the threshold has well-pronounced minima are 
found to be in very good agreement with numerical find- 
ings, see Eqs. (|32| ) and (f33|). Additionally, the analysis 
predicts that the threshold amplitude remains finite as 
the dissipative constant is vanishing, which is also con- 
firmed by the numerical results. 
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